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[i] The impact of climate change between 2000 and 2095 SRES A2 climates on surface 
ozone (0) 3 and on 0 3 source-receptor (S-R) relationships is quantified using three coupled 
climate-chemistry models (CCMs). The CCMs exhibit considerable variability in the 
spatial extent and location of surface 0 3 increases that occur within parts of high NO x 
emission source regions (up to 6 ppbv in the annual average and up to 14 ppbv in the season of 
maximum 0 3 ). In these source regions, all three CCMs show a positive relationship between 
surface 0 3 change and temperature change. Sensitivity simulations show that a combination 
of three individual chemical processes — ( i) enhanced PAN decomposition, (ii) higher water 
vapor concentrations, and (iii) enhanced isoprene emission — largely reproduces the global 
spatial pattern of annual-mean surface 0 3 response due to climate change (A 2 = 0.52). 

Changes in climate are found to exert a stronger control on the annual-mean surface 0 3 
response through changes in climate-sensitive 0 3 chemistry than through changes in transport 
as evaluated from idealized CO-like tracer concentrations. All three CCMs exhibit a similar 
spatial pattern of annual-mean surface 0 3 change to 20% regional 0 3 precursor emission 
reductions under future climate compared to the same emission reductions applied under 
present-day climate. The surface 0 3 response to emission reductions is larger over the source 
region and smaller downwind in the future than under present-day conditions. All three 
CCMs show areas within Europe where regional emission reductions larger than 20% 
are required to compensate climate change impacts on annual-mean surface 0 3 . 

Citation: Doherty, R. M., et al. (2013), Impacts of climate change on surface ozone and intercontinental ozone pollution: 
A multi-model study, J. Geophys. Res. Atmos., 118, doi: 10. 1002/jgrd. 50266. 


1. Introduction 

[ 2 ] Changes in climate are expected to influence future 
levels of surface ozone (0 3 ), a strong oxidant which has 
adverse impacts on health and ecosystems. Surface 0 3 is a 
local and regional pollutant which typically has peak 
episodes in spring and summer due to photochemical pro- 
duction. The long lifetime of 0 3 with respect to 


intercontinental transport timescales means that its influence 
on air quality can also be considered global [Akimoto, 2003; 
Holloway et al., 2003], and rising background 0 3 levels 
recorded at long-term measurement stations [e.g., Parrish 
et al., 2012] are of concern to policy makers. To determine 
the effectiveness of 0 3 precursor emission controls, the 
impacts and uncertainty in the 0 3 response to climate change 
need to be evaluated. 
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[3] Changes in climate at global and regional scales will 
modify the chemical environment and pollutant lifetimes 
and hence the concentrations of pollutants over source 
regions and over downwind continents [Task Force on 
Hemispheric Transport of Air Pollution (TF-HTAP), 
2011], Changes in climate may also affect meteorological 
transport processes and hence alter the export and import 
of pollution. Previous studies suggest that the response of 
O3 to climate change in polluted regions differs from that 
in remote regions [Murazaki and Hess, 2006]. In polluted re- 
gions, a positive 03-temperature relationship has been 
reported from both observational and model studies [Jacob 
and Winner, 2009]. Bloomer et al. [2009] define a climate 
penalty factor as the slope of the 0 3 -temperature relationship 
for percentiles of hourly O3 measurements over 2 1 years, 
which ranges between 2.2 and 3.2 ppbv/°C across the rural 
eastern US depending on emissions of NO x (NO+N0 2 ). 
Coupled climate-chemistry model (CCM) studies also show 
increased surface O3 associated with climate change in parts 
of major emission regions [e.g., Hauglustaine et al, 2005; 
Murazaki and Hess, 2006; Racherla and Adams, 2008; Wu 
et al., 2008; Royal Society, 2008; Jacob and Winner, 
2009]. Wu et al. [2008] define a climate penalty more broadly 
as the need for stronger emission controls to achieve a given 
air quality standard in a future climate. They find a climate 
penalty (O3 increase due to climate change) of 2-5 ppbv in 
summer daily maximum 8 hour O3 in the northeastern United 
States for climate change in the 2050s compared to the 2000s 
following the SRES A IB scenario. They attribute this 
increased O3 to a number of meteorological variables includ- 
ing temperature (a rise of 1-3°C) and mid-latitude cyclone 
frequency as well as increased biogenic isoprene emission. 

[4] The influence of climate change on O3 and its precur- 
sors occurs through multiple processes [e.g., Jacob and 
Winner, 2009; Isaksen et al., 2009; Fiore et al., 2012], 
Changes in temperature and water vapor alter the 
chemical environment and therefore affect the rates of chem- 
ical reactions that create and remove O3. Many chemical 
reaction rates increase with temperature, e.g., methane and 
non-methane hydrocarbon (NMVOC) oxidation rates, and 
lead to increased O3 production. In particular, thermal 
decomposition of peroxyacetylnitrate (PAN), a major reser- 
voir species for long range transport of the 0 3 precursors 
NO x and HO x (OH + H0 2 ), increases strongly with increas- 
ing temperature. Hence, increases in temperature will 
decrease the lifetime of PAN and contribute to reduced 
export of NO v (total oxidized nitrogen) and thus alter the 
long-range transport of 0 3 pollution [e.g., Schultz et al., 
1 998]. Increased water vapor in a future warmer atmosphere 
will lead to increased O3 destruction and shorter O3 
lifetimes [Johnson et al., 1999]. This is a robust feature of 
a previous multi-model study [Stevenson et al., 2006]. 
Amongst other influences on the tropospheric budget of 
0 3 , this may cause a reduction in the contribution of Asian 
emissions to background O3 over the United States 
[Murazaki and Hess, 2006; Lin et al., 2008]. However, in 
highly polluted regions, increased water vapor has compet- 
ing effects on O3 production [Jacob and Winner, 2009], 

[5] Climate change may also modify the future chemical 
environment through changes in natural emissions. Isoprene 
is a major 0 3 precursor under high NO x conditions [Trainer 
et al., 1991 \ Jacob and Winner, 2009]. Biogenic emission of 


isoprene increases strongly with temperature [Guenther 
et al., 1995, 2012], CCMs that employ interactive tempera- 
ture-sensitive emission schemes simulate net O3 production 
increases with higher isoprene emissions in high NO x 
regions but decreases in net O3 production in low NO x envi- 
ronments [e.g., Racherla and Adams, 2008; Zeng et al., 
2008]. Racherla and Adams [2008] find enhanced isoprene 
emission in future climate to be the dominant cause of 
increased summer O3 chemical production in the eastern 
United States. However, Ito et al. [2009] and Fiore et al. 
[2012] highlight that the sign of the 0 3 response to temper- 
ature and climate change depends on the assumption of the 
amount of recycling of NO x from isoprene nitrates. Further- 
more, the extent to which C0 2 inhibition of isoprene 
emissions in a future higher C0 2 climate may offset temper- 
ature-driven emission increases is unclear but may be 
substantial [Rosenstiel et al., 2003; Arneth et al., 2007; 
Heald et al., 2009; Young et al., 2009]. Other factors such 
as drought and changes in land cover and land use will mod- 
ify the spatial pattern and magnitude of isoprene emission 
[Sanderson et al., 2003; Wu et al., 2012, Guenther et al., 
2012], Another key uncertainty is how dry deposition 
changes in a warmer climate [Andersson and Engardt, 
2010; Wu et al., 2012] and in particular in regions where 
wanning is likely to be accompanied by drying, such as in 
the sub-tropics [Held and Soden, 2006]. Lastly, changes in 
cloud extent and properties as well as in precipitation will 
affect O3 production from lightning and also influence 
photolysis rates and wet deposition of nitric acid (HN0 3 ), 
the main NO x sink. 

[6] Meteorological transport pathways for pollutants may 
also be modified in a future climate through changes in 
synoptic and convective transport. A number of studies have 
suggested future decreases in synoptic-scale circulation 
frequency leading to increased summertime surface O3 pol- 
lution episodes over the eastern United States and Europe 
[Mickley et al., 2004; Forkel and Knoche, 2006; Murazaki 
and Hess, 2006; Leibensperger et al., 2008; Wu et al., 
2008], These changes in transport generally favor reduced 
export of pollutants from source regions. However, Racherla 
and Adams [2008] also find increasing 0 3 episodes, but 
neither these authors nor Lang and Waugh [2011] find 
evidence of future changes in synoptic-scale circulation. 

[7] This paper focuses on quantifying the impacts of 
climate change and their associated uncertainty on surface 
0 3 over its source regions and over downwind continents, 
hereafter termed source-receptor (S-R) relationships as in 
Fiore et al. [2011, 2009] and [TF-HTAP, 2007, 2011], 
It builds on the approach used in multi-model experiments 
coordinated by the Task Force on Hemispheric Transport 
of Air Pollutants (TF-HTAP). The Task Force, set up to 
inform the Convention on Long-range Transboundary Air 
Pollution (CLRTAP), designed model experiments to quan- 
tify source-receptor relationships by performing emission 
perturbations over four world regions: North America, 
Europe, East Asia, and South Asia [Fiore et al., 2009; 
TF-HTAP, 2007, 201 1], In this study, the model experiments 
are performed within the context of climate change between 
the 2000s and 2100s as simulated under the SRES A2 green- 
house gas emissions scenario that yields global surface 
temperature changes of ~3K between the two periods. Here 
three CCMs, evaluated in the previous TF-HTAP studies 
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for 0 3 and its precursors [ Sanderson et al., 2008; Shindell et 
al., 2008; Fiore et al., 2009; Jonson et al., 2010], are used to 
quantify the impact of climate change on surface O3 and on 
its S-R relationships. This study examines the effect of cli- 
mate change isolated from the effect of changing future 
emissions that are addressed in Wild et al. [2012]. Hence, 
the same (present-day) anthropogenic emissions are used in 
the simulations for “present-day” and “future” climates. 

[8] The methodology is described in section 2 and the ef- 
fects of climate change on surface O3 and its precursors are 
discussed in section 3. Section 4 analyzes 03-temperature - 
NO x relationships. Individual key chemical processes 
influenced by climate change are then examined as well as 
the influence of climate change on transport (sections 5 
and 6). S-R relationships for O3 are the focus of section 7, 
and the influence of individual chemical processes on O3 
S-R relationships are outlined in section 8. Section 9 com- 
pares the impact of climate change to the impact of anthro- 
pogenic emission reductions of 0 3 precursors on surface 
0 3 . Overall findings are then presented in section 10. 

2. Methodology 

[9] Model integrations were performed with three CCMs 
to examine (a) the influence of climate change on surface O3 
and its precursors and (b) the influence of climate change on 
intercontinental transport and O3 S-R relationships for four 
major emission regions. An additional set of experiments 
was performed with one CCM to quantify the relative impor- 
tance of different chemical processes influenced by climate 
change that impact surface O3 and its S-R relationships. 

2.1. Models and Climate and Emission Reduction 
Simulations 

[10] The influence of climate change on surface O3 and O3 
S-R relationships is simulated by three CCMs which comprise 
chemistry transport models coupled to atmospheric general 
circulation models (AGCMs): GISS-PUCCINI-ModelE (here- 
after GISS-PUCCINI), STOC-HadAM3, and UM-CAM. The 
three CCMs have been widely used [e.g., Fiore et al., 2009]. 


Their grid resolution and chemical, transport, and deposition 
schemes are outlined in Table SI in the auxiliary material. 

[11] To simulate present-day and future climate, the 
AGCMs were driven by sea-surface temperatures (SSTs) and 
sea-ice distributions from previous coupled ocean-atmosphere 
model integrations that were forced by greenhouse-gas emis- 
sions from the SRES A2 emission scenario [Nakicenovic 
et al., 2000]. The SRES A2 scenario was chosen as it lies at 
the upper end of the magnitude of future projected greenhouse 
gas emission trends and therefore represents a relatively large 
climate change signal with concomitant effects on chemistry 
and transport. In 2100, the SRES A2 scenario has a radiative 
forcing relative to the pre-industrial (1750) of about 8.1 W 
m 2 [IPCC, 2001], similar to the latest IPCC Representative 
Concentration Pathways (RCP) 8.5 scenario [Meinshausen 
et al., 2011], which has a radiative forcing of 8.5 W m 2 . 

[12] The GISS-Model E and HadCM3 GCMs were used to 
provide SSTs and sea ice for the GISS-PUCC1N1 and 
STOC-HadAM3/UM-CAM CCMs, respectively. To ensure 
that future changes in 0 3 and O3 S-R relationships can be at- 
tributed to climate change rather than inter-annual variability, 
5 years of simulations were performed. Two base-case integra- 
tions were carried out using meteorology from the AGCMs for 
2001-2005 and 2095-2099 (Table 1). The 5 year annual-aver- 
age increase in global-mean surface temperature and specific 
humidity between present-day and future periods averaged 
across the three CCMs was 3. OK (range 2.8-3.4K) and 19% 
( 1 8-2 1 %), respectively. These values are used for perturbation 
simulations in section 2.3. For GISS-PUCCINI, the 5 year av- 
erage surface temperature and humidity changes between 
2000-2005 and 2095-2099 were compared to differences in 
5 year averages separated by 100 years sampled from a 2000 
year long unforced control run to provide a comprehensive 
measure of internal model variability in the absence of 
climate change. Over much of the world, the surface tempera- 
ture and humidity changes from the climate change runs ex- 
amined here are 5 and 30 times greater, respectively, than 
those due to the unforced internal model variability (Figure 
SI). It is concluded that the changes in surface temperature 
and humidity between the 5 year simulations analyzed here 


Table 1. Experiments Performed and Their Description Including Number of Simulations and Models Performing Simulations 


Experiment Name Used in This Study 

HTAP Experiment 
Name 

Description 

Number of Models 

2000base 

FCl 

2000-2005 climate based on GCM SSTs 

All 3 models: GISS- 
PUCCINI, STOC- 
HadAM3, UM-CAM 

2095base 

FC2 

2095-2099 climate based on GCM SSTs 

All three models 

2000em_NA, EU, EA, SA 

FC3NA, EU, EA, SA 

Regional emission reductions of NO x , CO, NMVOCs 
for 2000-2005 climate for four world regions 

All three models 

2095em_NA, EU, EA, SA 

FC4NA, EU, EA, SA 

Regional emission reductions of NO x , CO, NMVOCs 
for 2095-2099 climate for four world regions 

All three models 

2000PAN a 


2000-2005 climate; added +3K to PAN decomposition 
rate: PAN + M = CH 3 C00 2 + N0 2 + M 

STOC-HadAM3 

2000H20 3 


2000-2005 climate; multiplied the water vapor 
concentration by 1 19%. 

STOC-HadAM3 

2000ISO a 


2000-2005 climate; added +3K to isoprene emission 
scheme 

STOC-HadAM3 

2000COM 3 


2000-2005 climate; added +3K to PAN decomposition 
rate and isoprene emission scheme and multiplied the 
water vapor concentration by 1 19%. 

STOC-HadAM3 


“For all these experiments, an additional simulation was performed with reduced emissions over the NA source region, i.e., 2000PANem_NA, etc. 
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are due to climate change and that for these simulations, 5 years 
is long enough to capture the climate change signal between 
2000 and 2095. 

[ 13 ] To isolate the effects of climate change on O3 and its 
S-R relationships, anthropogenic emissions and methane 
concentrations were held fixed at 2001 values for both 
present-day and future simulations, as in earlier TF-HTAP 
experiments [Fiore et al., 2009]. For GISS-PUCCINI, the 
methane concentrations were fixed at the surface at 1760 
ppbv, while for STOC-HADAM3 and UM-CAM, methane 
concentrations were fixed at 1760 ppbv throughout the 
model domain. Global emissions totals for NO x , CO, and 
NMVOC used in these simulations that include total, anthro- 
pogenic, biomass burning, and natural categories are given 
in Table S2, and the spatial distributions are shown in Figure 
S2. The three CCMs used similar anthropogenic and 
biomass burning emissions. Anthropogenic emissions for 
the year 2001 were largely based on the EDGAR3.2 dataset 
[Olivier and Berdowski, 2001], and biomass burning emis- 
sions were based on van der Werf et al. [2003] (Table S2; 
also for a comparison with RCP 2000 emissions). 

[ 14 ] Although anthropogenic and biomass burning emis- 
sions are unchanged between present and future, all the 
models use interactive lightning NO x emissions according 
to Price and Rind [1992, 1994] and Price et al. [1997] 
which are sensitive to changes in cloud top height in a future 
climate. The CCMs simulate future increases in lightning 
NO,- emissions ranging from 0.75 to 1.49 Tg N/yr (11-26%), 
and these increases are most prominent over S. America, 
Africa, and S. Asia (Figure S3). The STOC-FladAM3 model 
employs an interactive climate-sensitive isoprene emission 
scheme [Guenther et al., 1995] and gives a 22% global-mean 
increase in isoprene emissions due to climate change. Monoter- 
pene emissions which may also be sensitive to climate were not 
interactive in the STOC-HadAM3 simulations. Isoprene 
emissions remained fixed in the other two CCMs, reflecting 
uncertainty in the extent to which CO 2 inhibition of isoprene 
emissions in a future higher CO 2 climate may counteract tem- 
perature-driven increases. For example, Heald et al. [2009] find 
that increases in future isoprene emissions due to a projected 
warmer climate are entirely offset by including the CO 2 inhibi- 
tion effects. In a future 2090s climate under the SRES A2 
climate scenario with corresponding atmospheric CO 2 levels, 
Young et al. [2009] find either small increases or decreases 
(depending on location) in isoprene emissions when CO 2 
inhibition effects are included. However, there may also be an 
indirect effect of higher CO 2 that enhances isoprene emission 
through increased leaf area index (LAI) [Guenther et al., 2012], 

[ 15 ] These experiments termed “2000base” and 
“2095base” have the same setup as the control SRI simula- 
tions reported in Fiore et al. [2009] and TF-HTAP [2011], 
apart from differing climate conditions (Table 1). An 
addition to these simulations was the inclusion of artificial 
CO-like tracers in STOC-HadAM3 and UM-CAM. These 
tracers are emitted from anthropogenic CO sources and have 
a first-order decay lifetime of 50 days [Shindell et al., 2008; 
Fang et al., 2011; Schultz et al., 2013, in preparation]. One 
CO tracer was emitted from each HTAP source region. 
These tracers enable diagnosis of how changes in transport 
from source regions affect the distributions of trace gas 
species with similar lifetimes (such as CO and O3) between 
present day and future. The differences between the 


2000base and 2095base experiments are used to gauge the 
potential impacts of climate change on surface O3 and its 
precursors as well as on tracer transport. 

[ 16 ] To study the impact of climate change on O3 source- 
receptor (S-R) relationships, a further set of simulations was 
performed by the three CCMs (Table 1). These simulations 
focused on intercontinental transport between four major 
emission regions: North America (NA), Europe (EU), East 
Asia (EA), and South Asia (SA) as defined in Fiore et al. 
[2009], (hereafter termed the HTAP emission source re- 
gions) depicted in Figure 2a. For each of the four source 
regions, simulations were performed in which anthropogenic 
emissions of the 0 3 precursors NO x , NMVOCs, and CO 
were simultaneously reduced by 20%. These simulations 
were performed for both the 2000s and 2095 climates 
and are labeled by the region where emissions were reduced 
(e.g., 2000em_NA). Differences between the 2000em and 
2000base simulations for each source region give an esti- 
mate of the response of 0 3 and its precursors to a 20% emis- 
sion reduction over that region in a 2000 climate. The 20% 
regional perturbation was chosen to produce a clear O3 re- 
duction but yet to allow near-linear scaling to other perturba- 
tion sizes up to about 60% [Fiore et al., 2009; Wild et al., 
2012], 

[ 17 ] Since methane concentrations were fixed in all exper- 
iments, it is the short-term response of O3 to the changes in 
emissions that is diagnosed rather than any longer term O3 
response due to changes in methane lifetime resulting from 
OH changes [e.g., West et al., 2009; Fiore et al., 2009]. 

2.2. Model Evaluation for Present-Day Surface O 3 

[is] Fiore et al. [2009] compared simulated 0 3 in year 
200 1 from 2 1 models to observations over the United States, 
Europe, and Japan. Here Figure 1 shows how present-day 
5 year mean (2000-2005) surface 0 3 simulated in the 
“2000base” experiment by the three CCMs used in this 
study (shown as solid colored lines) compares with observa- 
tions and with the results from the 2 1 models in Fiore et al. 
[2009]. The three CCMs used here were also included in Fiore 
et al. [2009] (shown as dashed colored lines in Figure 1), but 
note that GISS-PUCCINI used different driving meteorology 
(NCEP reanalysis in Fiore et al. [2009] versus meteorology 
from a free-running GCM in this study) so small differences 
may be expected. Nonetheless, it can be seen that the 1 year 
(2001) and 5 year (2000-2005) mean 0 3 results from these 
three CCMs are fairly similar (Figure 1). 

[ 19 ] Overall, the seasonality of O3 simulated by the CCMs 
for the different locations lies well within the range 
simulated by the full set of models. The results from GISS- 
PUCCINI are nearer the lower end of the simulated O3 range 
in some regions, while STOC-HadAM3 on some occasions 
overestimates O3 compared to observations. GISS-PUCCINI 
exhibits a low 0 3 bias particularly in summer months for 
the Mediterranean, the Central and South Europe, and the 
SW United States, although there is an improvement in the 
5 year mean compared to the 2001 results for SW United 
States. However, its lower summer values in NE and SE 
United States are more in line with observations. The 
NMVOC emissions used by GISS-PUCCINI are lower than 
for the other two CCMs used here (Table S2) and compared to 
the ensemble-mean value of 630 Tg C yr -1 from Fiore et al. 
[2009]. This may in part explain the low O3 simulated by 
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Month of 2001 

Figure 1. Monthly average surface O3 observations (ppbv) for the year 2001 (solid black line; verti- 
cal bars represent one standard deviation) reproduced from Fiore et al. [2009, Figure 2], which de- 
scribes full details of the observations (CASTNET, EMEP, and EANET). The multi-model ensemble 
mean of monthly average O3 for 2001 from 21 models is shown with a black dashed line, and results from 
individual models are shown with gray lines. Monthly average surface O3 from the three CCMs applied in 
this study (GISS-PUCCINI: blue, STOC-HadAM3: red, and UM-CAM: green) are shown for 2001 as in 
Fiore et al. [2009] (dashed colored lines) and for the 5 year mean from the 2000base experiment (solid col- 
ored lines). 


GISS-PUCCINI. UM-CAM typically simulates 0 3 to within 
one standard deviation of the observations except in central 
Europe (Figure lc) and for summer and autumn in SW and 
W United States. STOC-HadAM3 typically overestimates 
O3 in summer in NE and SE United States and in the Great 
Lakes compared to observations, like the multi-model ensem- 
ble mean. Overall, the current three CCMs provide a represen- 
tative sample of the full set of models used in Fiore et al. 
[2009] in tenns of their present-day O3 simulation in 200 1 as 
shown in Figure 1. It is noted that the transport in STOC- 
HadAM3 and UM-CAM is similar since these two CCMs 
use the same AGCMs (section 2.1) but that the chemistry 
schemes are entirely independent (Table SI). 

2.3. Chemistry-Focused Sensitivity Simulations 

[20] To aid interpretation of the overall O3 response to 
climate change which results from a number of competing ef- 
fects, sensitivity experiments were perfonned with STOC- 
HadAM3 to isolate the 0 3 response to key chemical reactions 
and processes particularly sensitive to climate. Three prominent 
effects of climate change on surface O3 summarized by Jacob 


and Winner [2009] in their comprehensive review were investi- 
gated: enhanced PAN thermal decomposition (simulation 
2000PAN), enhanced water vapor concentrations (2000H20), 
and enhanced isoprene emissions (2000ISO) along with the 
combined effect of all three (2000COM); see Table 1. 

[21] These 5 year sensitivity simulations were carried out rel- 
ative to the 2000base simulation. The difference between these 
simulations and 2000base provides an estimate of the change in 
surface 0 3 due to each individual process. The similarity be- 
tween the spatial pattern of change in surface O3 to each 
perturbed process (and the combination of processes) and the 
surface 0 3 response due to climate change is quantified using 
a Pearson correlation (since a linear relationship was found 
between the two spatial datasets in all cases). The intention is 
to assess the relative importance of each process in different 
regions and the extent to which the combination of these three 
effects reproduces the net climate change effect on surface O3 

[22] In the “2000PAN” experiment, the average change in 
global-mean surface temperature of +3K across the three 
CCMs due to climate change (as derived from multi-model 
mean 2095base-2000base) was added throughout the model 
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domain to the temperature used in the calculation of the 
PAN decomposition rate: 

PAN + M^CH 3 C00 2 +N0 2 + M (Rl) 

[ 23 ] The decomposition rate increases very rapidly with 
temperature. For the reaction rate coefficients utilized in 
STOC-HadAM3, an increase in surface temperature from 
287K to 290K decreases the lifetime of PAN from ~4 to 

2.5 h, and at 425 hPa, an increase in temperature from 

2 5 OK to 25 3 K decreases the PAN lifetime from ~6 to 3 
months. The 2000PAN and the other sensitivity simulations 
do not account for the temperature and humidity change 
varying with altitude. However, the zonal-mean temperature 
and humidity increases due to climate change in the CCMs 
lie within the range of 3-5K and 15-30%, respectively, 
and are largely uniform with altitude in the lower and middle 
troposphere. The temperature increases due to climate 
change are about IK larger over land than ocean (Figure 
S4), and hence, using the global-average temperature change 
over the model domain in these sensitivity simulations will 
underestimate the effect of climate change due to tempera- 
ture over land and overestimate the effect over the oceans. 
Humidity increases due to climate change are more uniform 
over land and ocean (Figure S4). 

[ 24 ] In the “2000H20” simulation, the global average 
+19% surface increase in specific humidity was imposed 
throughout the troposphere. Increasing water vapor has a 
strong effect on O3 destruction through the reaction: 

0('D) + H 2 0 — > 20H (R2) 

[ 25 ] Increasing HO x also increases O3 destruction by 
direct reaction with O3. R2 is the primary source of OH rad- 
icals, and an increase in OH can either enhance O3 formation 
(through oxidation of NMVOCs, CH 4 , and CO) or suppress 
it (through increased N0 2 to nitric acid (HNO3) conversion) 
[Jacob and Winner, 2009]. 

[ 26 ] In the “2000ISO” experiment, a +3K temperature 
change was applied to the calculation of isoprene emissions 
from vegetation (the vegetation distribution remains fixed 
between present and future). This simple approach does 
not account for the effects on future isoprene emissions of 
C0 2 inhibition of emissions, changes in land cover, soil 
moisture or cloud cover, or future changes in land-use pat- 
terns. In STOC-HadAM3, isoprene emissions are related to 
temperature based on Guenther et al. [1995]. The emissions 
also depend on photosynthetically available radiation 
(PAR), but the climate sensitivity of this is not investigated 
here. In the scheme used, isoprene emission increases rap- 
idly with temperature up to a maximum at 314K. A 3K in- 
crease in surface temperature from present day values pro- 
duces a 17% increase in global isoprene emission. 

[ 27 ] The individual perturbations (+3K to PAN decom- 
position and isoprene emission and +19% increase to 
water vapor) were applied together in the “2000COM” ex- 
periment to assess the linearity of the surface O 3 response 
to these processes acting simultaneously. Finally, a further 
set of simulations were performed using the above simula- 
tions as a baseline and applying 20 % anthropogenic 
emission reductions of O 3 precursors in the NA source re- 
gion only. 


3. Climate Change Impact on Surface O 3 
and Its Precursors 

[ 28 ] Annual-mean surface O3 distributions for present-day 
(2000-2005) in all three CCMs highlight the northern mid- 
latitude emission and outflow regions where surface O3 con- 
centrations are high (Figure 2, top row). However, the spatial 
patterns of surface O3 across the globe are somewhat vari- 
able across the three CCMs; e.g., only STOC-HadAM3 ex- 
hibits very high annual-mean surface O3 in the southern 
hemisphere, which may be partly due to the higher NMVOC 
natural emissions in this CCM. However, differences in 
model chemistry and transport are also important since nat- 
ural NMVOC emissions in UM-CAM are only about 25% 
lower than in STOC-HadAM3 (Table S2). In contrast, natu- 
ral NMVOC emissions in GISS-PUCC1NI are substantially 
lower than that in the other two models which may explain 
the lower surface 0 3 concentrations in GISS-PUCCINI over 
the northern mid-latitude continents. Lower surface O3 over 
the N. Atlantic in STOC-HadAM3 is likely due to differ- 
ences in transport and deposition over oceans. The high 
surface O3 concentrations over the Tibetan plateau and over 
Greenland in GISS-PUCCINI (also seen in UM-CAM) are 
most probably due to higher stratosphere-troposphere 
exchange (STE) than in STOC-HadAM3. Overall, while dif- 
ferences in emissions account for some of the differences in 
simulated O3 between the three CCMs, differences in the rep- 
resentation of model chemistry and transport processes are 
likely to be the main source of these differences. 

[ 29 ] The annual-average surface O 3 response to the warm 
climate in 2095-2099 relative to 2000-2005 shows character- 
istic features described in previous studies [e.g., Hauglustaine 
et al., 2005; Murazaki and Hess, 2006; Zeng et al., 2008]: 
reduced surface O3 concentrations in less polluted regions and 
enhanced surface O3 concentrations in some polluted areas in 
all three CCMs (Figure 2, bottom row). Very similar spatial 
patterns of surface 0 3 response (and patterns of statistical signif- 
icance) were found when the STOC-HadAM3 and UM-CAM 
simulations were extended to cover two 10 year periods 
(2000-2009 and 2090-2099) (Figure S5). These results further 
confirm that 5 years is sufficient in these simulations to capture 
the climate change signal in surface 0 3 in all the CCMs. 

[ 30 ] Higher water vapor concentrations lead to reduced 
surface O3 in less polluted regions, causing a consistent de- 
crease in background O3 over most of the Earth’s surface in 
the CCMs [Thompson et al., 1989; Johnson et al., 1999, 
2001; Murazaki and Hess, 2006], In polluted regions (delin- 
eated very approximately by the 500 ppt NO x contours for 
the 2000 climate in Figure 2), the response is more mixed. 
Increased surface O3 (up to 6 ppbv) occurs over considerable 
parts of the major emission source regions in GISS-PUC- 
CINI and STOC-HadAM3 (Figures 2d-2e) but is confined 
to small areas in the UM-CAM simulations (Figure 2f). In 
the northern mid-latitudes, southern Europe and northeastern 
United States are the only regions where all CCMs consis- 
tently simulate O3 increases. Overall, the three CCMs 
exhibit some areas of O3 increase in the high NO x regions, 
but the spatial patterns of O3 increase in the three CCMs 
exhibit considerable variability in terms of their location 
and extent. A multi-model study using regional chemistry 
transport models also shows substantial variability in climate 
change-induced patterns of surface 0 3 increase over the 
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Figure 2. Top row: 5 year annual-average surface 0 3 concentrations (ppbv) (2000base) for (a) GISS- 
PUCCINI, (b) STOC-HadAM3, and (c) UM-CAM. Bottom row: The difference in the 5 year mean 
annual-average surface O3 concentrations (2095base-2000base) for the same three CCMs (d— f). The 
500 pptv contour of NO x surface concentrations for the 2000base simulation (thick black line) is used as an 
approximate indicator of polluted regions. Hatched areas, denoted by the + symbol, indicate where results are 
not significant at the 0.05 level as evaluated with a Student t-test using 5 years of data for the 2095 and 2000 
climate simulations. The HTAP source regions are depicted in Figure 2a: NA (15°N-55°N; 60°W-125°W), 
EU (25°N-65°N; 10°W-50°E), EA (15°N-50°N; 95°E-160°E), and SA (5°N-35°N; 50°E-95°E). 


United States [ Weaver et al., 2009]. Dilution of emissions 
over polluted regions due to the relatively coarse model 
resolution may result in a lesser sensitivity to climate change 
than would be given by higher resolution models. The 
strong O3 increase over the S. Hemisphere continents in 
STOC-HadAM3 is likely primarily due to increased iso- 
prene emission in these regions. Increased lightning NO x 
emissions over these regions may also influence surface O3 
concentrations especially for GISS-PUCCINI (Figure S3). 
The strong O3 increase simulated by GISS-PUCCINI over 
Tibetan plateau is likely related to enhanced STE in the 
future climate (Figure 2d). It is noted that the areas of O3 
increases from southern hemisphere continental outflow 
regions depicted for GISS-PUCCINI are not significant. 

[31] The season of maximum surface O3 for the 2000 
climate generally agrees with the season of maximum 
surface O3 derived from 2001 results (T. Nagashima, per- 
sonal communication, 2012), although in some northern 
parts of the mid-latitude continents, GISS-PUCCINI shows 
maximum surface 0 3 in winter due to the influx of strato- 
spheric air. The spatial patterns of the season of maximum 
surface O3 for the present-day within the four HTAP 
emission source regions vary typically between spring and 
summer across the CCMs; only in parts of the United States 
and in southern EU is summer (JJA) consistently the season 
of maximum O3 (Figure S6). The three CCMs again all 
consistently show O3 decreases during the season of 


maximum O3 but with areas of O3 increase in polluted north- 
ern mid-latitude regions. Increases in surface O3 due to cli- 
mate change in the season of maximum surface O3 reach 
up to 14 ppbv in parts of the HTAP source regions (Figure 
S6). Like the annual-mean surface O3 response, in the sea- 
son of maximum surface O3, there is considerable variability 
in the spatial patterns of 0 3 increase between the three 
CCMs due to climate change. Again, parts of the United 
States and southern Europe show consistent 0 3 increases 
across the three CCMs. In a multi-model regional study, 
Langner et al. [2012] also find consistent mean surface O3 
increases over April-September in southern Europe. 

[32] Overall, additional future emission controls would be 
needed to achieve a targeted level of O3 concentrations in 
many areas within polluted regions, e.g., southern Europe 
and northeastern United States, since annual-average O3 
concentrations are up to 6 ppbv or 10% higher (up to 14ppbv 
higher for the season of maximum surface O3) in 2095 than 
in 2000 due to climate change alone. However, the consider- 
able spatial variability and variation in the magnitude of 
simulated O3 increases in these CCMs preclude the quantifica- 
tion of broad regional O3 abatement targets. Generally, larger 
O3 increases of greater areal extent are simulated by STOC- 
HadAM3 (Figures 2 and S6); which is the only model to 
include a temperature-sensitive isoprene emission scheme 
(section 2.1). Langner et al. [2012] also find that out of four 
regional chemistry transport models, the model with the 
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Table 2. Regional Average Change in 5 Year Annual-Average Surface Meteorology and Chemistry Variables (2095base-2000base) 


Surface Variable 

Future (2095) Minus Present-Day (2000) Change in 

Annual-Mean Variable Averaged Over a HTAP Region 

NA region 

EU region 

EA region 

SA region 

Temperature 

4.4 (3.8-5.4)°C 

4.7 (4.1-5.4)°C 

4.2 (3.7-4.7)°C 

4.6 (3.8-5.0)°C 

Specific humidity 

19.2 (19.1-19.2)% 

20.6 (20.5-20.6)% 

23.9 (20.7-27.0)% 

24.6 (19.3-29.8)% 

Precipitation 

-4 (-10 to +0.6)% 

-17 (-8 to -22)% 

12 (7-14)% 

22 (20-25)% 

Fraction land area with 0 3 increases 

30 (2-47)% 

33 (10-46)% 

26 (2-53)% 

25 (3-67)% 

o 3 

-2.2 (-1.4 to -3.4)% 

-1.0 (-2.0 to +0.4)% 

-1.6 (-2.9 to +0.2)% 

-3.0 (-8.2 to +2.2)% 

PAN 

-30 (-24 to -42)% 

-32 (-22 to -44)% 

-25 (-20 to -34)% 

-37 (-35 to -39)% 

OH 

12 (2-18)% 

11 (7-15)% 

13 (9-16)% 

9 (5-12)% 

HN0 3 

0.8 (-5.6 to +5.0)% 

2.7 (1. 0-5.2)% 

8.0 (-3.1 to 19.3)% 

7.9 (0.2 to 22.2)% 

NO x 

-1.8 (-0.1 to -4.3)% 

-6.0 (-3.7 to -8.1)% 

-0.1 (-2.7 to 1.2)% 

-0.7 (—2.4 to +2.0)% 

HNO3 wet deposition 

0.4 (-9 to +7)% 

-17 (-14 to -19)% 

14 (11-22)% 

20 (12-28)% 

C 5 H 8 emission 3 

40% 

31% 

51% 

10% 

Net 0 3 production (P-L) 

13 (4-19)% 

19 (7-33)% 

7 (0.4-16)% 

3 (0.3-6)% 

0 3 dry deposition 

— 3.4(— 0.2 to -9.4)% 

-0.05 (-1.2 to +1.1)% 

-1.1 (-4.8 to +2.2)% 

0.1 (-10.3 to +11.0 )% 


Only land grid boxes are included in the averaging. Values represent the mean value across the three CCMs, with the individual model ranges given in 
parenthesis. 

“This variable is only available for STOC-HadAM3. 


largest isoprene emission change has the greatest sensitivity in 
its O3 response to climate change. Enhanced isoprene emis- 
sions in high NO x emission regions (where O3 production is 
VOC limited) can substantially augment O3 levels [e.g., 
Racherla and Adams, 2008] and is likely to be the dominant 
mechanism producing future O3 increases in the STOC- 
HadAM3 CCM. This is investigated further in section 5. 

[33] The influence of climate change on the O3 chemical 
environment was also investigated for the four HTAP emis- 
sion source regions (Table 2). Regional annual-average 
values were calculated within each HTAP region for land 
only. In each HTAP region, the fraction of land showing sur- 
face O3 increases is typically around one third, averaged 
across the three CCMs (Table 2). Overall, regional-wide 
annual-average surface O3 generally decreases slightly 
(1-3%) (Table 2) in all four HTAP regions. 

[34] Sillman and Samson [1995] suggest that 0 3 increases 
with temperature in polluted regions largely because of an 
increase in PAN decomposition (Rl). The production of 
PAN ties up NO x and reduces the source of peroxy radicals, 
which then can be subjected to long-range transport in the 
form of PAN [Murazaki and Hess, 2006]. Across the four 
regions, higher temperatures in the future (4.2-4.7°C) yield 
substantial decreases in PAN (25-37%) (Table 2) (in agree- 
ment with Hauglustaine et al. [2005]) and increases in 
hydroxyl radical (OH) concentrations. Higher water vapor 
mixing ratios (19-25%) also elevate OH (R2) (9-13%), 
which promotes O3 formation but also converts NO2 to 
HNO3 that is rapidly rained out to remove NO x and suppress 
further O3 formation [Jacob and Winner, 2009]. Accord- 
ingly, HN0 3 generally increases (1-8%), leading to a reduc- 
tion in the lifetime of NO x . Hence, unlike some previous 
studies [e.g., Murazaki and Hess, 2006], here PAN de- 
creases are typically accompanied by small NO x decreases 
or little change (0. 1-6%) (Table 2) rather than NO x increases 
across these source regions. This finding is however consis- 
tent with Racherla and Adams [2008] who report NO x de- 
creases (and HNO3 increases) associated with lower NO : 
NO2 ratios in high 0 3 regions across the eastern half of 
the United States in a future climate. Across the four HTAP re- 
gions, lower NO : NO2 ratios (0.5^4%) are simulated in future 
both in STOC-HadAM3 and UM-CAM (the only CCMs this 


ratio is available for). Changes in wet deposition of HNO3 vary 
in sign depending on region and are not consistently related to 
changes in HNO3 concentrations across the four regions but 
typically reflect the changing patterns of rainfall (Table 2). 

[35] In the warmer future climate, increases in isoprene 
emission in the four regions range from 10 to 50% in 
STOC-HadAM3, further promoting 0 3 formation in these 
high NO x emission regions. This result is reported in numer- 
ous other studies [e.g., Racherla and Adams, 2008; Zeng 
et al., 2008; Jacob and Winner, 2009]. Despite surface O3 
decreases, in all three CCMs, the changes in 0 3 chemistry 
described above lead to enhanced 0 3 chemical production 
and loss in the source regions such that net chemical produc- 
tion of O3 (chemical production (P)-chemical loss (L)) in- 
creases by 3-19% (Table 2). Another potentially important 
0 3 budget term is 0 3 dry deposition. However, changes in 
O3 deposition in these four source regions vary in sign across 
the CCMs, probably reflecting differences in model deposition 
schemes (e.g., inclusion of stomatal conductance) and near- 
surface meteorology (Table SI). The regional surface 0 3 re- 
sponse to climate change therefore appears dominated by the 
lowering of background 0 3 concentrations in a future climate 
which outweighs the increased regional net O3 chemical pro- 
duction (P-L). The lower regional-average surface O3 and 
higher O3 loss rate also imply a decrease in the O3 lifetime 
near the surface. Racherla and Adams [2008] also find in- 
creases in surface net chemical production and decreases in 
surface O3 lifetime over the eastern United States. 

[ 36 ] In summary, temperature and water vapor changes 
perturb NO v , HO x , and O3 chemistry through a number of 
complex interactions described above. The relative contribu- 
tions of temperature-dependent PAN decomposition (Rl) 
and isoprene emission changes, as well as water vapor (R2) in- 
creases to surface O3 change, are assessed in section 5. 

4. The Relationship Between Surface 0 3 
and Surface Temperature Change 

[37] In this section, the variability in surface O3 increase 
projected by the different CCMs due to climate change is 
explored in further detail. High-03 events in high emission 
regions often show a strong relationship with temperature 
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due to either chemistry processes or O3 precursor emissions 
sensitive to climate (section 5) or associated with stagnation 
episodes [Jacob and Winner, 2009; Rasmussen et al., 2012]. 
Here, the dependence of the relationship between the change 
in surface O3 with temperature and underlying NO x concen- 
trations is quantified for the three CCMs for the three mid- 
latitude HTAP regions: NA, EU, and EA. The analysis 
covers the months May to September which are typically 
the months of peak O3 for these regions (Figure S6). 

[ 38 ] Due to the substantial sub-regional variability in O3 
change over the large HTAP regions (Table 2), these three 
regions were divided into four sub-regions with equal 
latitude and longitude ranges. The SA region was not included 
because of its differing seasonal cycle of surface 0 3 (e.g., min- 
ima in JJA). Sub-regional averages for present-day NO x 
(2000-2005), A0 3 (2095-2099 minus 2000-2005), and AT 


GISS-PUCCINI 



STOC-HadAM3 



UM-CAM 



(2095-2099 minus 2000-2005) were calculated from all the 
grid cells over land within the appropriate sub-region. The 
relationship between AO3 and AT across all sub-regions, 
months, and years was then examined (sample size = 5 months 
x 5 years x 12 sub-regions = 300) (Figure 3). Positive rela- 
tionships were found for all three CCMs, with the climate pen- 
alty (the increase in surface O3 per 1°C increase in tempera- 
ture) ranging between 0.34 and 1.20 ppbv (Figure 3, left 
panels), although the correlation coefficient only reached 
r=0.45. The slopes were lowest for GISS-PUCCINI which 
also exhibits considerably lower mean NO x levels than the 
other CCMs and highest for STOC-HadAM3. 

[39] To assess the influence of NO x on AO3 /AT, the data 
were further divided into percentile ranges based on NO x 
concentrations in each CCM. The relationship between 
AO3 and AT was examined for the <25th, 25th-50th, 


GISS-PUCCINI 





Figure 3. Left-hand panels: Regional-average relationship between monthly AO3 versus AT (2095-2000 
climates) for each CCM for the three major mid-latitude source regions (NA, EU, and EA) for the peak 0 3 
months of May-September over all 5 years in each climate period (black line). Each point represents 
1 month in one of 5 years in each sub-region The data are color-coded by percentile category based on 
NO x concentration in the 2000s (orange = NO x concentration in the 50th-75th percentile; bright red = NO x 
concentration above the 75th percentile; dark red = NO x above the 95th percentile). The linear fit for A0 3 
versus AT was calculated for data that fell into each percentile, and each curve uses the same color coding. 
The slope value using all the data points is given at the top right of each panel, and the regional average 
NO x concentration (pptv) is given at the bottom right. Squares show the >95th percentile values that 
emanate from the EU region. Right-hand panels: the change in the slope of AO3/AT and correlation coeffi- 
cient (r) for each NO x percentile 0-25th, 25-50th, 50th-75th, 75th-100th, and 95th-100th with the x axis 
showing the NO x concentrations mid-way across the respective percentile range. Note the different x axis 
ranges for each CCM. 
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50th-75th, >75th, and >95th percentiles. In all three CCMs, 
all NO x values in the >95th percentile category, bar one 
point, are from the EU region (Figure 3, left panels). 

[40] No clear relationship between the gradient of AO3 and 
AT, and NO x level is exhibited for GISS-PUCCINI. For the 
other two CCMs, the gradient of AO3/AT steepens progres- 
sively for data in the 50th-75th and higher percentile catego- 
ries. Accordingly, the correlation coefficient is greater for 
these higher percentile categories (Figure 3, right panels). 
The NO x concentrations in these higher percentile ranges 
are also substantially larger in STOC-HadAM3 and UM- 
CAM compared to GISS-PUCCINI. For the >95th percen- 
tile data, the slope of AO3/AT and the mean NO x 
concentration are much higher in STOC-HadAM3 than in 
UM-CAM although the small sample size at this highest 
NO x level precludes definitive conclusions. Largely similar 
results are found when the analysis is repeated with de- 
seasonalized data. 

[41] To test the robustness of these results, a set of linear 
mixed statistical models [Pinheiro and Bates, 2000] were 
constructed to examine the association between simulated 
A0 3 and AT from each of the three CCMs. The effects of 
NO x , month, and region were considered, as were potential 
interactions between these variables and AT. Statistical 
models that included or excluded these variables were 
compared using Akaike’s Infonnation Criterion (AIC) 
[ Burnham and Anderson, 2002]. All of the statistical models 
also contained random effects for year, location, a year-by- 
location interaction, and a year-by-month interaction in order 
to avoid pseudo-replication by accounting for spatial and 
temporal dependence [Hurlbert, 1984]. 

[42] AIC provides a way of assessing the performance of 
statistical models by balancing the goodness-of-fit of the 
model (as measured by the maximum value of log-likeli- 
hood function, l) against the complexity of the model (as 
measured by the number of parameters in the model, p ) 
and is defined to be AIC = —21 + 2p. Statistical models with 
lower AIC values have better performance than models with 
high AIC values; the statistical model with lowest AIC can 
be regarded as being the best supported model, but models 
with AIC values close to this (typically AIC values within 
two units of the best model) [Burnham and Anderson, 
2002] can also be regarded as having good empirical sup- 
port. The best supported models (the statistical models with 
the lowest AIC values for each CCM) are listed in Table 3. 

[43] For all three CCMs there is strong evidence of a rela- 
tionship between AO3 and AT. Alternate statistical models 
that excluded AT were very poorly supported according to 
AIC. Hence, variations in AO3 cannot be explained solely 
by differences between months, regions, and NO x levels. 
For GISS-PUCCINI, the optimal statistical model for AO3 
is based on AT as a predictor alone (Table 3). For STOC- 
HadAM3, the best statistical model indicates that the magni- 
tude of the AO3 /AT relationship is dependent upon the NO x 


Table 3. Best Statistical Model for AO3 Based On AIC Criteria 


CCM 

Best Statistical Model 

GISS-PUCCINI 

STOC-HadAM3 

UM-CAM 

ao 3 =a t 

AO 3 = AT x NO x class + AT x month + region 
A0 3 = AT X NO x class + region 


level and upon month (since the model contained interaction 
terms between AT and NO x percentile and between AT and 
month). For UM-CAM, the best statistical model also shows 
an interaction between AT and NO x percentile. The best 
statistical models for STOC-HadAM3 and UM-CAM also 
contain an effect for “region,” suggesting that overall AO3 
values differ between regions but do not include an interac- 
tion between region and AT (suggesting that the strength of 
the AT-AO3 relationship does not vary between regions). Al- 
ternative models with good empirical support for A0 3 from 
STOC-HadAM3 and UM-CAM all contained AT and a 
dependence on NO x percentile. For GISS-PUCCINI, no al- 
ternate model was well supported. In summary, the results 
suggest that the slope of A0 3 /AT is influenced by NO x level 
for UM-CAM (most strongly) and for STOC-HadAM3 but 
not for GISS-PUCCINI, which has substantially lower mean 
NO x levels. 

[44] Overall, the results from STOC-HadAM3 and UM- 
CAM suggest that the sensitivity of surface O3 increases to 
temperature change is highest in the EU region due primarily 
to its high NO x levels. However, the slopes derived from the 
three CCMs here are generally lower than those derived 
from perturbation studies over polluted regions (>2 ppbv 
per 1°C increase) [Jacob and Winner, 2009], but these typi- 
cally use higher order O3 metrics than monthly mean O3 and 
consider smaller-scale urban regions. Rasmussen et al. 
[2012] also find modeled slope values of monthly average 
daily O3/T within two out of three regions in the eastern 
United States are lower in summer months compared to 
observations (their Figure 6). Steeper O3/T and AO3/AT gra- 
dients with higher NO x levels have also been previously 
reported. Bloomer et al. [2009] find a larger climate penalty 
(in this case, the slope of hourly O3 /T within four regions in 
the eastern U.S.) prior to 2002 when NO x emissions were 
considerably larger than that in later years. Wu et al. 
[2008] report a decrease in the climate penalty, defined as 
AO3 between 2050 and 2000, over the United States (their 
Figure 9), under a reduced (-40%) NO x emission scenario. 
Therefore, the modification of AO3/AT by NO x percentile 
range, shown by two of the three CCMs, is in broad agree- 
ment with the findings of these above studies. 

5. Chemical Processes Determining the Climate 
Change Impact on Surface 0 3 

[45] Results from sensitivity simulations with the STOC- 
HadAM3 CCM were used to determine the relative impor- 
tance of different chemical processes, influenced by climate 
change, on surface O3 changes. These sensitivity perturba- 
tions applied the difference in global-mean surface tempera- 
ture (3K) or water vapor (19%) between the 2095 and 2000 
climates as spatially uniform changes (section 2.3). 

[46] A 3K increase in temperature applied to the PAN 
decomposition rate (Rl) alone (2000PAN simulation) 
produces contrasting surface 0 3 responses over land and 
ocean. Over land, surface 0 3 concentrations increase by 0-2 
ppbv in northern mid-latitude regions and by up to 6 ppbv in 
the tropics (Figure 4a). Changes in O3 precursor concentra- 
tions and 0 3 production for these simulations are shown in 
Table 4 (for comparison with Table 2). Enhanced PAN de- 
composition leads to substantial decreases in surface PAN 
concentrations across the HTAP regions (25-32%), slight 
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a) PAN decomposition 


b) H 2 0 increase 
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Figure 4. The difference in 5 year annual-average surface O3 concentrations (ppbv) relative to 2000base 
for the following perturbations: (a) +3K to PAN decomposition rate (2000PAN), (b) 19% increase in wa- 
ter vapor mixing ratio (2000H20), (c) +3K to isoprene emissions scheme (2000ISO), and (d) the com- 
bined effect (2000COM). Results are from the STOC-HadAM3 CCM only. 


increases in surface OH, and increases in HNO3 (2 — 4%). 
These regional-average changes are similar to those due to cli- 
mate change except that the OH response is much smaller. 
Hence, the PAN and HNO3 responses due to climate change 
can largely be attributed to increased PAN decomposition. Re- 
gional-average surface NO x decreases slightly (except in EA) 
as seen under climate change. Regional-average net chemical 
0 3 production (1-6%) and surface 0 3 concentrations (~1%) 
increase over land (Table 4). 


[47] Less PAN transported from the emission regions 
leads to lower 0 3 production downwind and to surface 0 3 
decreases of up to 1 ppbv over the tropical oceans 
(Figure 4a). Hence, the effect of temperature increases on 
PAN decomposition contributes both to surface 0 3 enhance- 
ment in some parts of the HTAP source regions (~20-50% 
of the 0 3 increases due to climate change) and 0 3 decreases 
over the remote oceans (Figures 4a and 2b). The variability 
in the global spatial pattern of the O3 response to climate 


Table 4. Regional Annual-Average Change in Surface O 3 , O 3 Precursors and Net O 3 Chemical Production Between the Perturbed and 
Base-Case Present-Day Simulations (2000PAN, H20, ISO Minus 2000base) 


Surface Variable 


Perturbed Minus Present-Day Change in Annual-Mean Variable Averaged Over a HTAP Region (%) 



NA Region 


EU Region 


EA Region 


SA Region 


PAN 

H20 

ISO 

PAN 

H20 

ISO 

PAN 

H20 

ISO 

PAN 

H20 

ISO 

o 3 

1.0 

-4.6 

2.8 

0.8 

-3.9 

3.2 

1.6 

-4.1 

2.8 

0.7 

-5.4 

1.3 

NO x 

- 0.2 

- 1.0 

- 1.8 

-0.4 

- 1.2 

- 1.8 

0.4 

- 1.1 

- 1.8 

-0.5 

-0.9 

0.3 

PAN 

-29 

- 0.8 

13 

-27 

-0.5 

12 

-25 

-0.4 

12 

-32 

-1.5 

8 

HN0 3 

3.2 

1.1 

- 2.0 

2.1 

0.9 

-0.5 

4.1 

1.0 

- 2.8 

2.1 

0.3 

- 0.8 

OH 

1.6 

5.8 

- 6.0 

1.1 

6.0 

-3.3 

2.1 

5.9 

-7.4 

0.9 

5.7 

- 2.6 

Net 0 3 Prod.(P-L) 

5.9 

2.8 

7.3 

2.9 

2.0 

8.3 

5.1 

1.0 

5.6 

1.0 

-3.6 

0.4 

CsHg emission 



20 



18 



31 



6 


Only land grid boxes are included in the averaging. 
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change explained by the enhanced PAN decomposition 
alone is appreciable: r 2 = 0.3. In a further sensitivity simula- 
tion, +3K was applied to all chemical reaction rates. In this 
case, the O3 response was very similar in pattern to that from 
PAN decomposition (Rl) alone, but slightly smaller, with 
O3 enhancements in northern mid-latitude source regions 
not typically exceeding 1 ppbv (Figure S7). This suggests 
that increased temperature augments both O3 production 
(by enhancing PAN decomposition and OH oxidation of 
CH 4 and NMVOCs) and loss (by direct reaction of 0 3 , e. 
g., O3+NO2). Overall, enhanced PAN decomposition (Rl) 
is the primary driver of the O3 response and its spatial pat- 
tern over land and ocean. 

[48] The increase in water vapor concentrations alone, in 
2000H20, gives rise to increased 0 3 destruction through 
(R2). Surface 0 3 reductions of 1-2 ppbv occur globally 
and reach 3 ppbv across tropical locations (Figure 4b). In 
source regions, this decrease in surface 0 3 partly offsets 
the O3 increases caused by increased PAN decomposition. 
Across the four HTAP regions, the 19% increase in water 
vapor concentrations leads to uniform OH increases of 
~6% (Table 4). Changes in net O3 chemical production are 
generally positive but range from —4 to +3% depending on 
region. Surface O3 concentrations decrease by 4-5% across 
these source regions. Therefore, higher water vapor concen- 
trations in the 2095 climate relative to the 2000 climate are 
likely to be the primary cause of both the increase in OH 
and the decrease in surface O3 concentrations at the 
regional level (Tables 2 and 4). Over ocean regions, the 
increase in water vapor explains ~50% of the surface 0 3 
decrease due to climate change. The variability in the global 
spatial distribution of O3 response to climate change 
explained by enhanced water vapor is similar to that due to 
enhanced PAN decomposition (r 2 = 0.29). 

[49] A 3K increase in temperature on isoprene emissions 
alone leads to O3 increases over both land and ocean with 
the largest increases in the HTAP regions (~2ppbv) 
(Figure 4c). Enhanced isoprene emission only produces sur- 
face O3 decreases in small areas of Amazonia and Central 
Africa where isoprene emissions are already high and NO x 
levels are low. Across the four HTAP regions, isoprene 
emission increases by 6-31% (Table 4). Isoprene, as a 
NMVOC, is a major O3 precursor under high NO x condi- 
tions [Jacob and Winner, 2009] and also a precursor for 
PAN. Thus, isoprene affects the partitioning among oxidized 
nitrogen species, shifting the balance from HN0 3 towards 
PAN [. Horowitz et al., 1998; Fiore et al., 2011]. Accord- 
ingly, across the four HTAP regions, enhanced isoprene 
emission substantially increases surface PAN concentrations 
(8-13%) and slightly reduces HN0 3 . Surface OH concentra- 
tions are also reduced (3-7%) as noted in other studies 
[Racherla and Adams, 2008; Zeng et al., 2008; Fiore 
et al., 2011], and this contributes to reduced HN0 3 
[Racherla and Adams, 2008]. There are increases in 
regional-average net chemical production of 0 3 (up to 8%) 
and surface 0 3 (1-3%). These surface 0 3 increases account 
for ~40% of that due to climate change. Therefore, in these 
HTAP regions, the effect of temperature on isoprene 
emissions makes a larger contribution to increased surface 
O3 than enhanced PAN decomposition (Table 4 and 
Figures 4a and 4c) for the STOC-HadAM3 CCM. However, 
these O3 increases are larger than those simulated by Fiore 


et al. [201 1] with a different model. Over the NA region, a 
+3K increase in surface temperature leads to a 20% increase 
in regional-average isoprene emission and a 2.8 ppbv 
increase in annual-average surface O3 in STOC-HadAM3 
(Table 4). Fiore et al. [2011] find that a uniform 20% iso- 
prene emission increase over the NA region only increases 
monthly surface 0 3 by up to 0.6 ppbv (see their Figure 4). 

[50] Over the remote oceans, increased isoprene emission 
also leads to enhanced 0 3 , due to higher PAN concentrations 
[Pfister et al., 2008; Zeng et al., 2008, Fiore et al., 2011], 
Under climate change, these O3 increases partly offset the 
O3 decreases arising from increased water vapor and PAN 
decomposition. The insensitivity of isoprene emissions to 
temperature in the GISS-PUCCINI and UM-CAM simula- 
tions may therefore explain the smaller magnitude and areal 
extent of 0 3 increases due to climate change in high NO x 
regions than in the STOC-HadAM3 simulations (section 
3). However, isoprene nitrate chemistry [Fiore et al., 2005; 
Horowitz et al., 2007; Ito et al., 2009; Archibald et al., 
2010] and isoprene emission responses to increased temper- 
ature and atmospheric CO2 concentrations, as discussed in 
section 2.1, remain highly uncertain. 

[51] These sensitivity simulations illustrate the different 
influences on background O3 (through PAN and H2O 
changes) as well as on local O3 levels (through OH and net 
O3 chemical production changes) that contribute to the over- 
all continental-scale O3 response. The full climate change 
effect on surface 0 3 is likely to be more complex than these 
three individual chemical effects. However, the combination 
of these effects on surface O3 (Figure 4d) yields similar 
spatial patterns of change to that produced due to climate 
change (Figure 2b), r 2 value of 0.52. Furthermore, the sum 
of the surface O3 response to the individual effects is very 
similar in spatial pattern and in magnitude (to within ± 0.2 
ppbv) to the surface 0 3 response to the combination of 
effects. However, the magnitude of the combined 0 3 
response is about 40% lower in the four HTAP regions than 
the response due to climate change (Figure 2d). One 
reason for this is that the regional-average surface tempera- 
ture (3.7-5.4K) and humidity (19-30%) increases over the 
HTAP regions are substantially higher than the global-mean 
surface changes that were applied here (+3K for temperature 
and 19% for water vapor). Further sensitivity simulations 
were performed using the maximum regional-mean changes 
from the HTAP regions (+5.4K for temperature and 30% for 
water vapor; see Table 2). The 0 3 responses were found to 
scale linearly from the global-mean change to the maximum 
regional-average changes to within ±0.2 ppbv for the three 
effects. Furthermore, the combined effect of the maximum 
regional-mean changes produces O3 increases of a similar 
magnitude to those due to climate change (Figure S8). 
Finally, transport processes will also influence background 
and regional 0 3 responses to climate change, and these are 
discussed in the following section. 


6. Transport Changes 

[52] To investigate the influence of differences in transport 
between present-day and future climates, passive CO tracer 
species were utilized (section 2.2). The CO tracer emitted from 
the NA region (Figure 5a) is the focus of the discussion here. 
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a) Multi-model mean surface CO tracer from NA b) Multi-model mean surface ACO tracer from NA 



Figure 5. (a) Model-mean 5 year annual-average surface CO tracer concentrations (ppbv) for a species 
with a 50 day lifetime emitted as for CO over the NA source region from the STOC-HadAM3 and the 
UM-CAM CCMs (ppbv) for 2000 climate, (b) The difference in model-mean 5 year annual-average sur- 
face CO-tracer concentrations (ppbv) between the 2095 and 2000 climates. Both CCMs use the HadAM3 
GCM driven by the same SSTs for their simulations. Dotted areas, denoted by the • symbol, indicate 
where results are significant at a 0.05 level as evaluated with a Student /-test using 5 years of data for 
the 2095 and 2000 climate simulations. 


[53] The STOC-HadAM3 and UM-CAM CCMs that 
include these CO-tracer species use the same climate model 
forced by the same SSTs for the present-day and future 
climate simulations (although there are different advection, 
convection, and boundary layer mixing schemes within the 
two chemistry transport models; Table SI). The model-mean 
annual-average response of the surface CO tracer from NA 
to climate change is displayed in Figure 5. In the 2095 
climate, there is less (typically ~l-5% and up to 10%) of 
the annual-mean CO tracer from NA remaining at the 
surface in much of the NA source region than in the 2000 
climate (Figure 5b), especially over the Great Plains region. 
This may suggest an enhancement of the low-level jet 
[Murazaki and Hess, 2006] and enhanced venting and ex- 
port from the boundary layer in the future in that area. How- 
ever, over the eastern and western United States and outflow 
regions, the CO tracer increases (Figure 5b) (0-10% and up 
to 25%) in the future suggesting reduced ventilation from the 
surface. Only the areas of CO-tracer increase are statistically 
significant at the 0.05 level (Figure 5) although at the 0.01 
level, the area of CO-tracer decrease is significant. Similar 
spatial patterns of surface CO-tracer response were found 
between 2000 and 2095 when the STOC-HadAM3 and UM- 
CAM simulations were extended to cover two 1 0 year periods 
(2000-2009 and 2090-2099) and the areas of statistical signif- 
icance at the 0.05 level included the areas of CO-tracer 
decrease in the Great Plains region. This suggests that the 
CO-tracer changes are robust and likely due to climate change. 

[54] Overall, there appear to be complex shifts in transport 
that affect the surface distribution of the annual-mean CO- 
tracer response to climate change over NA and its outflow 
region, rather than any large-scale spatial response. Similar 
results were also found for CO-tracer species emitted from 
the other three HTAP regions (Figure S9), with distinct pat- 
terns of adjacent areas of lower and higher surface CO tracer 
concentrations that suggest a shift in circulation within all 
the regions that extends across the regional boundaries 
between present day and future. It is unlikely that these 
shifts in transport patterns have a major role in influencing 
the spatial patterns of at least the annual-mean O3 response 
due to climate change. However, changes in transport may 
well be important when considering changes in higher 


moments of surface O3 such as changes above the 95% per- 
centile of 0 3 or in daily 8 h maximum surface 0 3 . 

[55] However, as these transport-related results are derived 
essentially from one GCM climate simulation, this precludes 
definitive statements about the robustness of the projected fu- 
ture changes in annual-mean CO tracer transport. Fang et al. 
[2011] consider the impact of climate change on CO-tracer 
transport using a tracer emitted globally with a 25 day life- 
time in the GFDL AM3 model. They find a vertical redistri- 
bution of CO tracer in the future climate, which is also 
found here. However, Fang et al. [2011] find that surface 
zonal-mean CO tracer concentrations increase in the tropics 
and mid-latitudes in the future relative to present day. How- 
ever, STOC-HadAM3 and UM-CAM simulate uniform 
increases in zonal-mean CO tracer only above -800 hPa 
with mainly decreases below this pressure level. Therefore, 
CO-tracer transport changes may vary and depend on the 
representation of shallow convection processes in the 
individual models. 

[56] In contrast to the mixed CO-tracer response to climate 
change over the HTAP regions found here, a number of pre- 
vious studies have highlighted reduced boundary layer 
venting in a future climate over parts of the United States 
[Mickley et al., 2004; Murazaki and Hess, 2006; 
Leibensperger et al., 2008; Wu et al., 2008] and W. Europe 
[Hauglustaine et al., 2005; Forkel and Knoche, 2006], 
Racherla and Adams [2008] suggest that different methodol- 
ogies for analysis of synoptic-scale circulation changes may 
produce different results. Further process-based studies 
using a range of CCMs and methodologies are needed 
before more robust statements can be made concerning the 
effect of climate change on tracer transport. 

7. Climate Change Impact on Ozone S-R 
Relationships 

[57] In this section, the analysis is extended to consider 
0 3 S-R relationships and how these are modified due to 
climate change. Following the same approach as previous 
HTAP studies [Fiore et al., 2009; TF-HTAP, 2011], 0 3 
S-R relationships were quantified by considering the O 3 
response to a 20 % reduction of 0 3 precursors emissions 
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Figure 6. Impact of simultaneous 20% reductions of the 0 3 precursor emissions (NO x , NMVOCs, and 
CO) in the NA source region on 5 year annual-mean surface 0 3 concentrations (ppbv) for 2000 climate for 
(a) the GISS-PUCCINI CCM, (b) STOC-HadAM3, and (c) UM-CAM (2000em_NA— 2000base). The 
HTAP regions are depicted in Figure 6a). 
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Figure 7. The reduction in regional-average annual-mean surface 0 3 in the receptor region from 20% 
reductions in 0 3 precursor emissions (NO x , NMVOCs, and CO) in the (a) NA, (b) EU, (c) EA, and (d) 
SA source regions for the 2000 (00) and 2095 (95) climates. MM01 = Multi-model mean and ranges with 
200 1 climate [Fiore et al., 2009] . ME00 = 3 model-mean and ranges for 2000 climate; ME95 = 3 model-mean 
and ranges for 2095 climate; for ME00 and ME95, results are 5 year mean averages. No land-sea mask 
is applied to the calculation of these area averages. 
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Figure 8. The difference in the impact of 20% reductions of O 3 anthropogenic precursor emissions over 
NA on 5 year annual-mean surface O3 (ppbv) for the 2095 climate compared to the 2000 climate for the (a) 
G1SS-PUCCINI, (b) STOC-HadAM3, and (c) UM-CAM CCMs (2095em_NA-2095base-2000em_NA- 
2000base). Hatched areas, denoted by the + symbol, indicate where results are not significant at the signif- 
icance level of 0.05 as evaluated with Student /-test using 5 years of data for the 2095 and 2000 climate 
simulations. Negative values represent a larger reduction whilst positive values imply a lesser reduction in 
the 0 3 response to emission reductions. 


of NO x , CO, and NMVOCs over a given source region, 
as compared to “no reduction” (section 2.1). The HTAP 
source regions discussed in the previous sections are now 
considered as both source and receptor regions for 0 3 and 
its precursor emissions (Figure 6 a). The NA source region 
is used for illustrative purposes in the next two sections. 

[58] For present-day climate, in response to a 20% 
emission reduction of O3 precursors over the NA region, 
all three CCMs simulate annual-mean surface O3 decreases 
of 0.5-2ppbv (~2 — 10%) in the source region and up to 
lppbv (~3%) immediately downwind (Figure 6 ). Similar 
surface O3 responses are found over the other three HTAP 
regions (Figures 7b-7d). The magnitude of the surface 0 3 
response simulated by GISS-PUCCIN1 is typically slightly 
less compared to the other two CCMs for all four HTAP 
source regions. The decreases in annual-mean surface O3 
over source and receptor regions simulated by all three 
CCMs for 2000 climate (ME00) (Figure 7) lie well within 
the multi-model range of results from the 1 5 models that par- 
ticipated in the sensitivity simulations for 2001 (MM01) 
(Figure 7) [Fiore et ah, 2009], 

[59] The difference in the annual-mean O3 response due to 
the 20% reduction of O3 precursors in the NA region 
between the 2000 and 2095 climates is depicted for the three 
CCMs in Figures 7a and 8 . When averaged over the large area 
of the HTAP regions, the difference in the surface O3 response 
to emission reductions between the two climate periods is 
small (Figure 7). However, the differences in the spatial pat- 
tern of surface O3 response to emission reductions between 
present-day and future climate is robust across the three CCMs 
(Figure 8 ) despite differences in their spatial patterns and mag- 
nitude of O3 response to climate change (Figures 2d-2f). 

[60] Across the three CCMs, the O 3 response to NA emis- 
sion reductions in the future climate is consistently larger 
over the source region and smaller downwind compared to 
present-day climate, both in terms of absolute and relative 
changes. The annual-average surface O3 response to the 
20% NA emission reduction is enhanced by up to 0.6 ppbv 
(2%) in the NA source region and reduced by up to 0.2 ppbv 
(0.3%) immediately downwind (Figure 8 ). The regionally 
averaged annual-mean surface O3 response to emission reduc- 
tions is larger by 0. 1 ppbv for the NA source region and 
smaller by -0.05 ppbv and -0.03 ppbv for the EU and EA 
receptor regions respectively (Figure 7a) in the future climate. 


[61 ] Equivalent emission reduction simulations for the 
three other source regions yield similar results, with a greater 
surface O3 response to a 20% emission reduction in the 
source region and a reduced response in downwind conti- 
nents in the future compared to present-day climate 
(Figures 7b-7d). Over the SA region, the effect of large-scale 
averaging over land and downwind ocean leads to cancelation 
effects; hence, there is little difference in the regional O 3 
response to emission reductions between the two climate 
periods (Figure 7d). In addition, the effect of emission reduc- 
tions in the SA source region has a smaller influence on its 
three receptor regions compared to the influence of the other 
source regions on their respective receptor continents 
(Figure 7). This reflects the smaller region size and the more 
southerly location of this region such that SA pollution 
remains isolated from mid-latitude air [Fiore et al., 2009]. 
Across all source and receptor regions, there is little inter- 
annual variation (< 0.02 ppbv in the source regions; not 
shown) in the difference in the surface O 3 response between 
the two 5 year climate periods. 

[62] The seasonality of the O3 response to the reduction of 
O3 precursors remains similar for the two climate periods in 
both the source and receptor regions (Figure S10) in the 
three CCMs. Over the source region, the EU shows the larg- 
est changes in magnitude of the O3 response between the 
two periods in all three CCMs. There is little evidence of a 
consistent change in the peak month(s) of intercontinental 
transport between source and receptor regions, or in the 
frequency of underlying transport (e.g., storm track and flow 
patterns) between months, or in the seasonality of chemical 
transformation processes (e.g., O3 lifetime) between the 
two climate periods. Thus, climate change has little influ- 
ence on the seasonality of the O3 response to 20% emission 
reductions. This result is consistent with Fiore et al. [201 1] 
who found the seasonality in the O3 response to 20% de- 
creases in NA anthropogenic emissions was unaltered when 
isoprene emissions were elevated. 


8. Chemical Processes Influencing the Climate 
Change Impact on Surface 0 3 S-R Relationships 
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[63] The climate-sensitive chemical processes discussed in 
section 5 are re-considered here in the context of how they 
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Figure 9 . The difference in the impact of 20% reductions of O3 precursor emissions over NA on 5 year 
annual-mean surface O3 (ppbv) relative to 2000base for the following perturbations: (a) +3K to PAN 
decomposition rate (e.g., 2000PANem_NA-2000PAN-2000em_NA-2000base), (b) 19% increase water 
vapor mixing ratio, (c) +3K to isoprene emissions scheme, and (d) the combined effect. Results from 
the STOC-HADAM3 CCM only. Negative values represent an enhanced reduction, whilst positive values 
imply a lesser reduction in the O3 response to emission reductions. 


influence the annual-average O3 response to O3 precursor 
emissions reductions over the NA region (Figure 9). 

[ 64 ] The difference in the surface O3 response to 20% NA 
emission reductions of O3 precursors between the sensitivity 
simulations with increased temperature applied to the PAN 
decomposition (Rl) alone (2000PANem_NA-2000PAN) 
and the base (2000em_NA-2000base) simulations (Figure 9a) 
shows a similar spatial pattern to that due to climate change 
(Figure 8b) (r 2 = 0.5). There is a greater 0 3 response to emis- 
sion reductions in the source region but a lesser response 
downwind (Figure 9a). Flowever, the difference in the 0 3 re- 
sponse is smaller than that due to climate change. Under 
higher water vapor concentrations alone, only a small part of 
the NA source region shows an enhanced surface O3 response, 
and there is a smaller O3 response elsewhere in the source re- 
gion and downwind (Figure 9b). This downwind 0 3 response 
to emission reductions for higher water vapor concentrations 
alone is ~50% of the magnitude of that due to climate change 
(Figures 8b and 9b). In contrast, under elevated isoprene 
emission arising from a uniform +3K increase in temperature, 
there is a greater O3 response to 20% anthropogenic emission 
reductions in both source and receptor regions (Figure 9c). 
These differences in the absolute O3 response largely reflect 
the changes in surface 0 3 concentration that occur under these 
different perturbation simulations. 

[ 65 ] The three processes thus have competing effects on 
the O3 S-R response to emission reductions in both the NA 
source and receptor regions. For the source region, the 
temperature effect on isoprene emission influences the surface 
O3 response more than PAN decomposition. In the downwind 


region, there is a reduced surface 0 3 response due largely to 
water vapor effects but also to increased PAN decomposition. 
The combined effect of the three processes produces a spatial 
pattern of O3 response to emission reductions which is similar 
to that due to climate change (Figures 9d and 8b) (r 2 = 0.76), 
but the magnitude is somewhat smaller. However, scaling to 
the maximum regional-mean temperature and humidity 
changes for these three effects yields O3 responses similar in 
magnitude to those due to climate change. 

9. The Impact of Climate Change Versus 
Emission Reductions on Surface 0 3 

[66] Given that climate change increases surface O3 in 
some areas within major emission source regions, a key 
question for policy makers is how effective current air qual- 
ity legislation will be in the future. The percentage change in 
region-wide 0 3 precursor emissions required to balance the 
localized impact of climate change (i.e., at each CCM grid 
point) on annual-average surface O3 concentrations over 
each HTAP region is shown in Figure 10. This is based on 
the parameterization of Wild et al. [2012] which provides 
an estimate of regional O3 responses based on regional pre- 
cursor emission changes scaled quadratically from the 20% 
emission reductions of NO x , CO, and NMVOCs applied in 
the HTAP studies. This approach reproduces the nonlinear 
behavior of regional O3 responses well for emission reduc- 
tions of up to about 60%, and therefore, the focus here is 
on changes up to this magnitude to minimize the associated 
uncertainty. For simplicity, we assume that the emissions of 
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Figure 10. The percentage reduction in region-wide 0 3 precursor emissions (NO x , NMVOCs, and CO) 
required to counteract the impact of climate change (2095base-2000base) on 5 year annual-average sur- 
face 0 3 concentrations over each HTAP region. This is calculated by dividing the impact of 2095-2000 
climate change on surface O3 by the effect of 20% regional O3 precursor emission reductions on surface 
0 3 and multiplying the result according to the parameterization of Wild et al. [2012] (e.g., by 20% for a 
linear scaling). Each row shows a different model, and each column shows a different source region. Neg- 
ative values (in blue) show where climate change leads to a reduction in O3 and therefore where emissions 
can increase to maintain present-day (2000s) O3 levels. Areas are shaded white when 20% regional 
emission reductions lead to increased 0 3 . 


each of the different O3 precursors are altered proportion- 
ally. Positive values (in red) show areas within the HTAP 
regions where regional emission reductions are necessary 
to overcome the increases in surface O3 due to climate 
change such that present-day O3 levels are retained. The 
spatial pattern of required regional emission change within 
each HTAP region varies considerably with CCM. This pre- 
dominately reflects the variation in the annual-mean O3 
response due to climate change rather than variability in 
the O3 response to regional 0 3 precursor emission reduc- 
tions. Although all HTAP regions show areas where regional 
emission reductions of more than 20% are required to bal- 
ance annual-average surface O3 increases due to climate 
change, this result is most consistent across the CCMs for 
the EU region. Within this region, there are large areas 
where regional emission reductions of more than 20% are 
required to balance the impact of climate change between 
2095 and 2000 on annual-average surface O3 concentrations 
(Figure 10), and reductions of more than 60% are required 
in some places. For the SA region, only the results from 
GISS-PUCCINI suggest that substantial regional emission 
reductions would be required to mitigate the localized 
effects of climate change across the Tibetan plateau (via 


enhanced STE) on the annual-mean 0 3 response (Figures 10 
and 2d). Larger regional emission reductions are generally 
required in STOC-HadAM3 to overcome climate change 
than in the other two CCMs. Again, this reflects the greater 
0 3 response to climate change in STOC-HadAM3 which 
may emanate from the temperature control on biogenic 
emissions that is only simulated in this CCM. 


10. Conclusions 

[67] This study examines the impact of climate change on 
surface O3 and its precursors and on intercontinental trans- 
port of O3 from major emission regions to downwind receptor 
continents. It further explores the relative contributions of 
key processes in producing this 0 3 response. Three CCMs 
(GISS-PUCCINI, STOC-HadAM3, and UM-CAM) simulate 
surface O3 concentrations for 2000 and 2095 climates 
as projected under the SRES A2 climate forcing scenario 
with current and reduced emissions of anthropogenic O3 
precursors. To isolate the effects of climate change, anthropo- 
genic and biomass burning emissions and methane concentra- 
tions are fixed at 2000 levels. 


17 



DOHERTY ET AL.: CLIMATE CHANGE IMPACTS ON 0 3 AND TRANSPORT 


[68] In the three CCMs, annual-mean surface O3 generally 
decreases due to climate change but increases in parts of the 
more polluted regions, consistent with previous studies. 
Surface 0 3 increases in high NO x regions are up to 6 ppbv 
in the annual average and up to 14 ppbv in the season of 
present-day maximum surface 0 3 . This climate penalty 
effect on surface 0 3 within the high NO x regions is highly 
variable across the CCMs in location and spatial extent. 
The largest areas and magnitudes of O3 increase are simu- 
lated by STOC-HadAM3, the only CCM that incorporates 
temperature-sensitive isoprene emission. Annual regional- 
mean surface 0 3 decreases slightly (1-3%) in conjunction 
with substantial PAN decreases ( — 30%) and elevated OH 
(~10%) and HNO3 levels across the four HTAP regions. 
This shift in NO v partitioning from PAN to HN0 3 under 
climate change is in agreement with previous studies [e.g., 
Zeng et al., 2008; Racherla and Adams, 2008], 

[69] The relationship between monthly (May-September) 
surface O3 change (AO3) and temperature change (AT) is 
quantified in order to explore how the climate penalty relates 
to NO x concentrations. Positive AO3/AT relationships are 
exhibited by all three CCMs and the gradient of AO3/AT 
steepens with increased NO x percentile class in UM-CAM 
(most strongly) and STOC-HadAM 3 . Consequently, in 
UM-CAM and STOC-HadAM 3 , the sensitivity of surface 
O3 increase to temperature change is greatest in the EU 
region, where all three CCMs simulate the highest NO x con- 
centrations. The lack of variation in AO3/AT with NO x per- 
centile class for GISS-PUCCINI may reflect its lower NO x 
concentrations. 

[70] The relative importance of three chemical processes 
that affect 0 3 chemistry and are sensitive to temperature 
and humidity changes (PAN decomposition, water vapor 
concentrations (hence HO x ), and isoprene emission) is ex- 
amined in further simulations using the STOC-HadAM3 
CCM. An increase in temperature acting only on the PAN 
thermal decomposition rate leads to surface O3 increases 
over land and decreases over ocean. Increased water vapor 
mixing ratios lead to O3 decreases everywhere. In contrast, 
enhanced isoprene emission due to higher temperature yields 
surface O3 increases over most locations. Given this poten- 
tial importance of future changes in isoprene emission fur- 
ther research is needed to quantify uncertainties in isoprene 
nitrate chemistry and the influence of temperature and CO2 
effects on isoprene emissions [e.g., Guenther et al., 2012], 
Other climate-sensitive emissions such as methane from 
wetlands [e.g., Shindell et al., 2006] and NO x from soils 
[e.g., Zeng et al., 2008] also merit further study. Although 
the climate change impact on surface O3 is likely to be more 
complex than the sum of individual chemical effects, it is 
found that the combination of these three effects largely 
reproduces the spatial patterns of annual-mean 0 3 response 
due to climate change ( R 2 = 0.52). 

[71] The impact of climate change on transport is evalu- 
ated with an idealized CO-like tracer species emitted from 
all four major source regions. Over all four HTAP regions, 
annual-mean CO-tracer concentrations exhibit distinct 
dipole patterns of increase and decrease in future compared 
to present-day climate that suggest subtle shifts in transport 
within the region and its area of outflow, rather than any 
consistent large-scale pattern. It is unlikely that these com- 
plex shifts in regional circulation exert a substantial control 


on the coherent spatial pattern of the annual-mean surface 
0 3 response to climate change; hence, the simulated 
annual-mean surface O3 response to climate change arises pre- 
dominately from changes in chemistry. However, the effect of 
climate change on transport may well depend on model 
representation of shallow convection processes. Further pro- 
cess-based studies are required to quantify how climate 
change affects the transport of chemical constituents. 

[72] The impact of climate change on intercontinental 
transport of 0 3 is also considered through source-receptor 
(S-R) relationships that quantify the impact of changing 
anthropogenic emissions in a source region on surface O3 
both within and beyond that region. A small but robust dif- 
ference in the annual-mean surface 0 3 response to a 20% 
reduction in 0 3 anthropogenic precursor emissions (NO x , 
CO, and NMVOCs) between the present day and future cli- 
mate is simulated by all three CCMs. The O3 response to re- 
gional emission reductions is larger over the source region 
and smaller downwind in the future climate over all four 
HTAP regions. There is no change in the seasonality of the 
O3 response to emission reductions in either the source or re- 
ceptor regions. The combination of the surface O3 response 
to emission reductions under (i) enhanced PAN decomposi- 
tion, (ii) higher water vapor mixing ratios, and (iii) elevated 
isoprene emission, relative to the O3 response to emission re- 
ductions under present-day climate, is similar to that due to 
climate change (R 2 = 0.76). 

[73] Overall, this study suggests that the dominant effects 
of mean climate change on annual-mean O3 and its intercon- 
tinental transport are via temperature and water vapor effects 
on the chemical environment rather than climate-related 
changes in transport. Changes in transport may well be im- 
portant when considering changes in peak 0 3 events. 

[74] The regional emission reduction required to balance 
the impact of climate change such that present-day annual- 
mean surface O3 levels are retained is calculated across the 
four HTAP regions. The required emission change within 
each HTAP region varies considerably with CCM; typically 
larger values are projected by STOC-HadAM3, the CCM 
with temperature-sensitive isoprene emission. All three 
CCMs consistently show areas within the EU region where 
substantial (>20%) regional emission reductions are 
required to balance climate change impacts on annual-mean 
surface 0 3 . This emphasizes that the impact of climate 
change must be accounted for in designing future emission 
policies at least for these four major emission regions. 
Langner et al. [2012] find that the spatial features of the 
change in daily maximum O3 produced by climate change 
are similar to those for summer time average O3. Hence, this 
current study focusing on annual and seasonal mean changes 
should be a useful guide for assessing the effects of climate 
change on policy-relevant O3 metrics. 
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